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Abstract 

Background: Saccharomyces cerevisiae strains isolated from natural settings form structured biofilm colonies that 
are equipped with intricate protective mechanisms. These wild strains are able to reprogram themselves with a 
certain frequency during cultivation in plentiful laboratory conditions. The resulting domesticated strains switch off 
certain protective mechanisms and form smooth colonies that resemble those of common laboratory strains. 

Results: Here, we show that domestication can be reversed when a domesticated strain is challenged by various 
adverse conditions; the resulting feral strain restores its ability to form structured biofilm colonies. Phenotypic, 
microscopic and transcriptomic analyses show that phenotypic transition is a complex process that affects various 
aspects of feral strain physiology; it leads to a phenotype that resembles the original wild strain in some aspects 
and the domesticated derivative in others. We specify the genetic determinants that are likely involved in the 
formation of a structured biofilm colonies. In addition to FLOll, these determinants include genes that affect the 
cell wall and membrane composition. We also identify changes occurring during phenotypic transitions that affect 
other properties of phenotypic strain-variants, such as resistance to the impact of environmental stress. Here we 
document the regulatory role of the histone deacetylase Hdalp in developing such a resistance. 

Conclusions: We provide detailed analysis of transcriptomic and phenotypic modulations of three related S. 
cerevisiae strains that arose by phenotypic switching under diverse environmental conditions. We identify changes 
specifically related to a strain's ability to create complex structured colonies; we also show that other changes, such 
as genome rearrangement(s), are unrelated to this ability. Finally, we identify the importance of histone deacetylase 
Hdalp in strain resistance to stresses. 
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Background 

In the natural environment, microorganisms do not exist in 
stable and optimal growth conditions. Instead, microbial 
populations are often challenged by harmftil external fac- 
tors. The formation of complex multicellular communities, 
such as biofilms or colonies, provides an advantage for sur- 
vival in the wild due to the cooperative behavior of cells 
and the establishment of common protective mechanisms 
[1]. Moreover, the adaptive responses of individual cells and 
their phenotypic heterogeneity, factors that enable efficient 
adaptation to a rapidly changing environment, are also im- 
portant characteristics that allow populations to thrive in 
natural settings. Phenotypic switching is a mechanism by 
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which phenotypic heterogeneity is achieved. This process 
generates cell phenotypes that are better adapted to a par- 
ticular environment, thereby allowing a population to react 
flexibly to environmental changes [2]. In yeast, phenotypic 
switching generates a diverse array of changes, leading to 
the emergence of colonies with altered morphologies [3]. 
The appearance of different colonies occurs more fre- 
quently than somatic mutations. Although relatively stable 
during subsequent passages, the switched variants retain 
the ability to revert to their original phenotype [4]. 

Phenotypic switching is a phenomenon that is commonly 
observed in pathogenic yeasts [5]. Switching occurs at 
infection sites to generate phenotypes that effectively evade 
the host immune response [6]. The new phenotypes also 
exhibit features that contribute to increased virulence [7,8]. 
A phenotypic switch is defined as a random event that 
is not necessarily induced by external factors; however. 
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environmental conditions can affect the frequency and 
direction of the switch. For example, certain strains of 
Candida albicans undergo a transition from a "white" 
to "opaque" colony phenotype and vice versa. Exposure to 
temperatures above 30°C increases the frequency of 
switching to the "white" phenotype [9]. Thus, white cells 
predominate in the blood stream at a physiological 
temperature of 37°C, whereas "opaque" cells more effi- 
ciently colonize the skin surface, an environment with a 
lower temperature [7]. Furthermore, anaerobic conditions 
and other environmental factors also mediate the transition 
to the "opaque" phenotype [10]. Thus, switching provides a 
mechanism of adaptation to certain niches and to a variety 
of physiological conditions. 

Wild Saccharomyces cerevisiae strains isolated from 
their natural habitat also exhibit phenotypic heterogeneity 
and the ability to switch between two or more different 
colony pheno types [11-15]. During laboratory cultivation 
on rich media, the switch is usually oriented toward the 
formation of less-structured colonies that differ in many 
features from their structured counterparts [14]. Such a 
switch, that we term "domestication" [13] results in the 
formation of strains that are stable during subsequent pas- 
sages on agar media and that form smooth colonies simi- 
lar to those formed by standard S, cerevisiae laboratory 
strains [13]. The formation of a domesticated derivative, 
BR-S, of the wild S, cerevisiae BR-F strain on a non- 
fermentable medium occurs efficiently with an average 
frequency of approximately 2-3% [13]. This frequency sug- 
gests that switching is caused by a regulated event rather 
than by random mutations. During cultivation of a BR-S 
strain under adverse conditions, we demonstrate the oc- 
currence of reverse phenotypic changes that lead to the 
re-appearance of cell clones that form more structured 
colonies. We termed these strains "feral" strains. By com- 
paring the transcrip tomes of the original wild BR-F strain, 
its domesticated BR-S derivative and the feral BR-RF strain 
(the strain that restored biofilm colony formation), we de- 
termine the genome-wide expression alterations involved 
in phenotypic changes. In addition, we specify the factors 
involved in the formation of structured biofilm colonies. 
Finally, we show that at least some individual strain pro- 
perties are under the control of epigenetic mechanisms 
and that recombination and genome rearrangements 
occur during phenotypic switching. 

Results 

Feral subclones derived from a BR-S strain under stress 
conditions form structured colonies 

To induce the conversion of a BR-S strain to a strain with a 
wild-type phenotype, the opposite process of domestication, 
we set up various stressful and long-term starvation condi- 
tions. We incubated the BR-S strain statically (i.e., without 
shaking) for several months in various media with a limited 



carbon source. During incubation in MM medium with 2% 
ethanol, the number of colony-forming units (CFU) in the 
suspension was monitored. In parallel, the morphology 
(structured versus smooth) of the arising microcolonies was 
determined on GMA plates (Figure lA). After inoculation, 
the static cell culture grew slowly until approximately the 
60th day, as indicated by the slowly increasing number of 
CFU. From approximately the 75th day, the number of CFU 
started to decrease, suggesting a gradual dying of part of the 
population (Figure lA). During this period, the number 
of cells forming smooth colonies decreased to a rate of 
about 5 X 10^ cells per ml of culture per day. Colonies with 
a structured morphology started to appear rarely among 
the smooth colonies on the 28th day or later. The fre- 
quency of their appearance reached approximately 0-7% of 
the CFU (Figure lA). Interestingly, between days 30 and 
110, the number of cells forming structured colonies in- 
creased to a rate of up to 400 to 2000 cells per ml of the 
culture per day in parallel cultures. Thus, at later stages of 
static cultivation (140 days), the proportion of structured 
colonies increased to as high as 26% (Figure lA). In some 
static cultivations the cells forming structured colonies 
reached as high as 40% of all CFU (not shown). Structured 
colonies also emerged in static cultivations with MM 
medium that contained lower or higher ethanol concentra- 
tions. In addition, colonies appeared in cultivations without 
any carbon source, in which ethanol was replaced with 
NaNOs. When ethanol was replaced with fructose, struc- 
tured colonies only appeared when the starting fructose 
concentration was 0.5% or lower; with a fructose concen- 
tration of 1 - 2%, no structured colonies were observed in 
cultivations lasting more than 200 days. 

The morphology of the emerging structured colonies 
was not uniform, varying from moderately structured to 
a fully structured colony phenotype. We classified the 
emerging colonies into four typical colony phenotype 
categories (Figure IB). Most of the isolated subclones 
were unstable; after re-plating, these colonies domesti- 
cated at a high frequency, thus preventing a more de- 
tailed characterization. The most stable subclones 
formed colonies of the IV^^ category, phenotypically re- 
sembling biofilm colonies of the original wild BR-F 
strain. One of these subclones, designated BR-RF, was 
used for further analyses. 

The BR-RF strain domesticates more efficiently than the 
original wild BR-F strain 

With the aim of monitoring the stability of the BR-RF 
strain, we compared the frequency of its domestication on 
fermentative and respiratory media (evaluated as the fre- 
quency of formation of smooth colonies) with that of the 
BR-F strain (Figure IC). In BR-F colonies, the percentage 
of domesticated subclones ranged from 0.05% (4-day-old 
colonies) to 2.6% (22-day-old colonies) when grown on 
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(See figure on previous page.) 

Figure 1 Origin and stability of the BR-RF strain. A. Upper panel, growtli curve of 3 parallel static cultures of the BR-S strain in MM-2% ethanol 
medium, expressed as CFU/ml. Lower panel, appearance of cells forming structured colonies (switched phenotype) expressed as a percentage of 
total colony number. A representative experiment of the three is shown. B. Colonies with structured morphology originating from the BR-S strain 
during cultivation in static culture. From the top, morphotypes I, II, III and IV (the latter named BR-RF). Scale bar= 1 mm. C. Domestication frequency 
of colonies of the BR-F and BR-RF strains. The upper graphs represent the average appearance of smooth colonies after the re-plating of BR-F (red columns) 
and BR-RF (blue columns) colonies relative to the total CFU of all re-plated colonies. The bottom graphs show the switching frequencies of individual 
colonies, BR-F (red columns) and BR-RF (blue columns). 



respiratory medium and from 0.6% (4-day-old colonies) to 
7.5% (22-day-old colonies) when grown on fermentable 
medium (Figure IC, upper graphs). In individual BR-F 
colonies, the frequency of domesticated subclones varied 
from 0% to 15% on respiratory and from 0% to 24% on fer- 
mentable medium (Figure IC, bottom graphs). In BR-RF 
colonies grown on respiratory GMA, domesticated sub- 
clones appeared with an average frequency of 2% (4-day-old 
colonies) to 21% (22-day-old colonies), with variation in in- 
dividual colonies from 0% to 37%. Similar to BR-F domesti- 
cation, the frequency of BR-RF domestication was even 
higher when the BR-RF colonies were grown on fermentable 



YEGA; domestication varied from 10% (4-day-old colonies) 
to 26% (22-day-old colonies), with variation in individual 
colonies from 2% to 53% (Figure IC). 

Morphology and physiology of colonies formed by the 
BR-RF strain 

The morphology of colonies (birds-eye view) formed by the 
BR-RF strain was identical to that of colonies formed by the 
original BR-F strain (Figure 2A). Two-photon excitation 
confocal microscopy (2P-CM) confirmed that the BR-RF 
colony was composed of an aerial part and of subsur- 
face roots attaching the colony to the substrate. This 
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Figure 2 Ultrastructure and physiological properties of colonies of the BR-RF strain and its ancestors. A. Colonies of S. cerevisiae BR-F, BR- 
RF and BR-S strains. Scale bar= 1 mm. B. Vertical transversal cross-section of an 3.5-day-old BR-RF colony stained with ConA-AF and visualized by 
2P-CM. Scale bar= 100 |jm. At the bottom, magnification of the central region of the BR-RF colony (left) and lower subsurface part of the colony 
(right). Scale bar = 10 pm. C. Physiological characteristics of BR-F, BR-RF and BR-S colonies. Left, amount of wet biomass per colony. Right, relative 
amount of dry mass. D. Localization of Flol 1 p-GFP and morphology of cells taken from 4-day-old BR-F and BR-RF colonies. Right, Nomarski contrast 
image of cells from 4-day-old BR-S colonies. The Flol 1 p-GFP production in these colonies was below the detection level. Scale bar = 5 pm. 
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composition is similar to the organization of the BR-F 
colony [16]. The colonies of both strains contained in- 
ternal cell-free cavities (Figure 2B). In contrast to BR-F 
colonies [16], BR-RF colonies were formed exclusively of 
oval cells, including the subsurface parts where chains of 
oval cells invaded the agar (Figure 2B) instead of the pseu- 
dohyphae that are typical of BR-F colonies. Thus, although 
the BR-RF strain gained the ability to form fully structured 
colonies, it did not revert its cell morphology and retained 
the oval cells typical of the BR-S strain (Figure 2D). 

Although they formed similarly structured colonies, BR-F 
and BR-RF strains may have differed in other physiological 
parameters. We therefore analyzed BR-RF colonies grown 
on GMA for prominent characteristics previously found to 
be different between colonies formed by the BR-F and BR-S 
strains [14]. The analysis revealed that in some aspects, 
BR-RF colonies behaved as an "intermediate" between 
BR-F and BR-S colonies. This intermediate behavior 
was mostly related to wet biomass accrual (Figure 2C) 



and colony size (Figure 2A), while the water content (rela- 
tive amount of dry mass) of BR-RF colonies was almost 
identical to BR-F colonies (Figure 2C). 

Another typical feature of BR-F colonies was the presence 
of high-molecular-weight glycosylated protein (HMWGP), a 
possible constituent of the extracellular matrix (ECM). 
The production of HMWGP was switched off in BR-S 
colonies [13]. BR-RF colonies restored production of 
this protein, but at a lower level than in BR-F colonies. 
In addition, SDS-PAGE showed that the mobility of this 
protein was shifted to lower MW. This observation sug- 
gested that the molecule was shortened or underwent less 
glycosylation (Figure 3A). The presence of HMWGP (an 
ECM constituent) and a high content of water, possibly 
trapped by ECM [14] implied that BR-RF colonies restored 
ECM production. To demonstrate presence of ECM we 
prepared a BR-RF-Pgali-GFP strain to monitor the ECM 
in the developing colony by 2P-CM, as previously de- 
scribed [16]. As in BR-F colonies [16], internal BR-RF 
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Figure 3 Visualization of tlie ECM, presence of HMWGP protein and FLOll mRNA in BR-RF colonies. A. Glycoproteins in extracts from BR- 
F, BR-S and BR-RF colonies stained with silver (Ag) and visualized by the ConA/peroxidase (ConA) method; arrows indicate HMWGP in BR-F (black 
arrows) and in BR-RF (red arrows) colonies. B. Visualization of the ECM in vertical cross-sections of BR-RF-pcAirCFP colonies using 2P-CM. Green, 
GFP fluorescence indicates areas in which the inducer (galactose) reached the cells; red, autofluorescence of all colony cells visible in areas where 
the ECM prevented the inducer from accessing the cells. Intact colonies were induced from the bottom by placing them for 5 h on agar soaked 
with 2% galactose [16]. C. FLOl 1 mRNA in 4- and 7-day-old BR-F, BR-RF and BR-S colonies. Right, stronger exposure revealed low FLOl 1 mRNA 
levels in BR-S colonies. The discontinuous line compares the length of the FLOll transcripts in BR-RF and BR-S cells. The RDN18 gene was used as 
an RNA loading control. 
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colony regions were also protected via a low-permeable 
ECM (Figure 3B). 

The production of the adhesin Flollp is indispensable 
for the formation of structured colonies [14,17], and 
Flollp production is switched off during wild strain do- 
mestication [13]. Both FLOll gene expression (Figure 3C) 
and Flollp protein production (Figure 2D) are recovered 
in the BR-RF strain. Interestingly, the FLOll mRNA of 
the BR-RF strain is longer than that of the BR-F strain 
(Figure 3C). 

Genome-wide expression differences between BR-F, BR- 
RF and BR-S colonies 

To obtain an overall view of the transcription character- 
istics typical of the biofilm colony phenotype, we isolated 
total RNA from 4-day-old and 7-day-old colonies 
formed by the BR-F, BR-S and BR-RF strains. We then 



compared the transcriptomes of the BR-RF versus BR-F 
strain, the BR-RF versus BR-S strain and the BR-F versus 
BR-S strain using a microarray (Additional file 1: Table SI). 
Mutual comparison of the transcribed genes enabled 
the identification of genes specifically activated or re- 
pressed in complex structured biofilm colonies (repre- 
sented by the BR-F and BR-RF strains) compared to 
smooth colonies (represented by the BR-S strain) 
(Table 1). However, the expression differences related 
to the structured biofilm colony phenotype were only a 
subset of the expression differences observed between 
the strains. Other extensive differences unrelated to the 
colony phenotype were observed between BR-RF and 
BR-F colonies, as well as between BR-RF and BR-S 
colonies. Some of these differences are summarized in 
Figure 4 and discussed below. This finding supports the 
data described above indicating that some properties of 



Table 1 Functional classification of genes up-regulated in particular colony phenotypes 





Structured biofilm colony morphotype 


Smooth colony morphotype 


Cell wall 


SVSl, ASP3-1r2r3r4, BSC1, 05^5 


SPS100, CTS1, F/.09, PESh PSAl, ZPSl 


Secretory 
transport & 
processing 


BSTl, LST7, PERI, CWH41, IV1NN9, YOL107W, 
APM4, GVP36, YIP4 


BTN2, NCE102, GYP7, m4, PEP12, SNC2, TIM17, TIM21 


Metabolism 


JLPl, MLSl, SNZl, SNZ2, TPIl — YMR226Q 
BNAl, ENo'h HPTl, MDH2, RHR2, SOLA 


DSFl, PIG2, STR3, CRS5, CUPl-l, EBP26, E0L2, ERE6, HXKl, MTDl, NESl, SEAl 


Amino ocids 


ARG3, ARG4, BATl, MET17, MET22, ARG1, 
GDH3, LAPS 


GCVl GCV2, GCV3, SHM2 


Lipids & fatty ocids 


INOl — ECU POTI, POXl, TESl, VPS66 


YNL144C, ERGS, ERG4, HFDl — AURl, CEMl, CSG2, SURl, SUR4 


Nucleotides 




RNR4 — ADEh ADE2, ADE8, ADEIS, PRSS 


Transport 


ITRl, PUT4, VHTl, PDR12, VBA3, AT02, 
CAN], GAPh MEP2, PIC2 


YFL054C HXT15, MUP3, PH084, THI74, ZRCl, HXT5, HXT17, PH089, YEL054C, 
HXTh HXT2, HXT3, HXT4, HXTG, HXTIO, , HXT12, HXT13, MALll, PH087, ZRT2 


Degradation 


HSPSS, MGRh NAS2, YSP3 


ATG8, BSD2, CICl, DEG16, GID8, JIDl, UMPl, RPN5, RPNIS, RPT2, SNAS, VID24, YPT53 


Stress response 


DDR48, GRX3 — TSAl, YBL036Q CTAl, S0D2, 
TSA2 


OLAl, HSP42, SPGl, SSA3, SSA4, APJ1, FES1, GPX1, HSP30, HSP78, SSE1, 

YNLWOW^ AAD6, CCT5, CPR6, GLOl, GSHl, HSP104, HSP150, MGEl, 0YE3, PSTl, 
RCN2, SEDl, STIl, TRX3, YDJIYGPI 


Regulation 






Translation & RNA 
processing 


BRR2, MTQ2, SLFl, DUS4, SPP381, IST3, 
MRS2, RPPIA, RPR2 


PUS2, SLHl, NGL3, RPA135, RSC30, YDR341Q CBE5, SR09, RPP2A, RPS30B, RPL35B, 
RRP3, RRP7, SRCl 


Signal transduction 


EARS 


GPA2, MTHl, MRKl, STDl, LST8, RGS2, YMR291W^ GLC8, HRKl, KICl, MTHl, STE7, 
YAKl 


Cell cycle & polarity 


SHE1 


PCL9, HBTl, TAHl 1, CLG1, CDC/3, CDC3, DSEh VHSl 


Transcription 


CUP9, SFG1,TEC1, YAPA 
DAL80, RDSl 


YFL052W, HMSI, RLM1, USV1, ZPR1, ALPHA2^ C//V5, GISh HAP4, OPIl, PH02, 
STP2, STP4 


- chromatin 
structure 


HHFl 


T0S8, I0C4, HHTh TMA2SJOD6 


General 
transcription 




IWRh RBA50, RPA1S5, SPT15, SUA7, TAEU, TEAl 


Transposons^ 


8 genes, 10 genes, 1 genes 




Unknown* 


7 genes, 5 genes, 17 genes 


8 genes, 4 genes, / / genes 


Others 




YLR164W, DIA1, DIA3, MRH1, IVIRPS17, IXRh PINS 



Genes only expressed differently on the 4th day are in standard letters, genes differently expressed on both the 4th and 7th day are in bold, and genes only 
expressed differently on the 7th day are italicized. 

*genes are listed in the Additional file 1: Table SI and discussed in the text. 
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Figure 4 Comparison of transcriptomes from cells forming BR-F, BR-RF and BR-S colonies. Venn diagrams comparing genes induced in 
particular strain comparisons as indicated in titles of A, B and C panels. Red circles, comparison of BR-F versus BR-S; green circles, comparison of BR-S 
versus BR-RF; and blue circles, comparison of BR-F and BR-RF (see also the legend at the bottom of the figure). 



BR-RF colonies differ from colonies formed by both BR-F 
and BR-S. 

Expression characteristics of complex biofilm colonies 

Comparison of the structured biofilm colony versus the 
smooth colony transcriptomes revealed 73 and 89 genes 
with significantly higher expression in 4- and 7-day-old 
structured colonies, respectively. In smooth colonies, the 



expression levels of 73 and 157 genes were increased in 
4- and 7-day-old colonies, respectively (Figure 4, Table 1, 
Additional file 1: Table SI). 

The structured biofilm colonies up-regulated genes in- 
volved in protein secretion and modification. Some of 
these genes, such as MMV9, CWH41, BSTl and PERI 
influence the structure and the composition of the cell 
wall and/or the membrane. In addition, the structured 
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colonies highly up-regulated genes for cell wall proteins 
such as BSCl, which is similar to the FLO 11, and ASPS 
genes. On the other hand, the structured colonies down- 
regulated the expression of other cell wall-related genes 
such as SPSlOO and CTSl as well as genes related to se- 
cretion and protein modification. These data indicate 
that structured and smooth colonies differ in the carbo- 
hydrate and protein composition of their cell wall 
(Table 1). 

The structured and smooth colonies also differed in 
their expression of metabolic genes and nutrient trans- 
porters. The structured colonies up-regulated amino 
acid metabolism genes and transporter genes, including 
genes involved in arginine and methionine metabolism. 
The structured colonies also expressed genes involved in 
phospholipid biosynthesis and fatty acid |3-oxidation, 
and they expressed the inositol transporter gene ITRl, 
In contrast, the few amino acid metabolism genes that 
were up-regulated in the smooth colonies were for the 
mitochondrial glycine decarboxylase complex. The 
smooth colonies also up-regulated a number of nucleo- 
tide metabolism genes, genes involved in copper and 
iron resistance and homeostasis and several genes re- 
lated to lipid metabolism. These included genes involved 
in ergosterol, sphingolipid and ceramide metabolism 
(ERGS, SURl, SUR4, CSG2), The transporter genes in- 
cluded those for the transport of hexoses and phosphate 
{PH084, PH087) (Table 1). 

In addition, the smooth colonies up-regulated a large 
group of genes related to the stress response (28 genes), 
including genes for various chaperones, genes involved 
in cell wall stress response/maintenance and certain 
genes involved in the oxidative stress response. The few 
stress response genes induced in structured colonies in- 
cluded oxidative stress response genes {GRX3, TSAl, 
TSA2, CTAl and S0D2) and genes involved in DNA 
damage response and repair {DDR48, YBL036C, DDIS 
and RADIO) (Table 1). 

The cells of structured and smooth colonies also dif- 
fered in their expression of a remarkable number of 
genes involved in signaling cascades (1 gene up- 
regulated in the structured and 13 genes in the smooth 
colonies), transcription regulators (6 genes up-regulated 
in the structured and 13 genes in the smooth colonies) 
and genes involved in translation and RNA processing (9 
genes up-regulated in the structured and 14 genes in the 
smooth colonies). The transcription factor genes up- 
regulated in the structured colonies included TECl and 
SFGl, which are involved in yeast filamentation and 
pseudohyphal growth. In contrast, the smooth colonies 
up-regulated RLMl, YFL052W and USVl, which are re- 
lated to cell wall integrity. They also up-regulated several 
genes related to chromatin structure and remodeling 
{T0S8, T0D6, TMA23, 10C4 and HHTl) (Table 1). 



In addition, the structured colonies predominantly 
expressed genes for transposons and a large group of 
genes with unknown function, some of which were lo- 
calized to the subtelomeric regions (Table 1). 

BR-RF colonies up-regulated expression of genes located 
at specific positions in chromosomes 

BR-RF colonies increased the expression of 244 genes 
and 379 genes in 4- and 7-day-old colonies, respectively, 
compared with BR-F colonies. BR-RF colonies reduced 
the expression of 259 genes and 509 genes in the 4- and 
7-day-old colonies, respectively, compared to BR-F col- 
onies. For the comparison between the BR-RF and BR-S 
strains, BR-RF colonies increased the expression of 276 
genes and 390 genes in 4- and 7-day-old colonies, re- 
spectively, compared to BR-S colonies (including typical 
"biofilm colony" genes, see above). BR-RF colonies re- 
pressed 220 and 428 genes in 4- and 7-day-old colonies, 
respectively, compared to BR-S colonies. In addition, a 
relatively large number of genes were up-regulated (103 
genes in 4-day-old and 150 genes in 7-day-old colonies) 
and down-regulated (58 genes in 4-day-old and 176 
genes in 7-day-old colonies) in BR-RF colonies when 
compared with both BR-F and BR-S colonies (Figure 4). 

"Positional Gene Enrichment" analysis revealed that 
many gene expression differences between the BR-RF 
strain and either the BR-F or BR-S strain (or both) in- 
volved the expression of genes localized to specific 
chromosomal regions (Figure 5) rather than genes ran- 
domly distributed in the genome. Most prominently, a 
high number of genes up-regulated in 4-day old colonies 
of BR-RF strain were localized to a 300-kbp region of 
the left arm of chromosome XV. This region included 
33% of all activated genes in the BR-RF versus BR-S 
strain, 32% of all activated genes in the BR-RF versus 
BR-F strain and 60% of all activated genes in the BR-RF 
versus both BR-S and BR-F strains. Smaller groups of 
genes up-regulated in the BR-RF strain versus either the 
BR-S or BR-F strains were localized to a 320-kbp region 
of chromosome IX , and a group of BR-RF strain genes 
that were down-regulated compared to the BR-F and 
BR-S strains was localized to a 240-kbp region of 
chromosome XII (Figure 5). In contrast, neither the gene 
expression differences between the BR-F and BR-S 
strains, nor the genes expressed specifically in structured 
biofilm colonies (i.e., those differently expressed between 
the BR-F and BR-RF strains versus the BR-S strain) were 
localized to specific chromosomal regions (not shown). 

As the higher level of gene expression related to dis- 
tinct chromosomal areas could be caused by DNA dupli- 
cations, we measured the copy numbers of 3 genes 
located in the left arm of chromosome XV (the genome 
region with higher expression in BR-RF), relative to 
those of 3 genes located in the right arm of the same 
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Figure 5 Position clustering of differently expressed genes on chromosomes XV, IX and VII. Positional gene enricliment analysis performed 
on the basis of expression comparison of 4-day-old colonies of BR-F versus BR-S, BR-RF versus BR-S and BR-RF versus BR-F. X axis, position of each gene along 
the particular chromosome, the position of the centromere is marked by a green circle on green bar on the bottom of each graph; Y axis, relative expression 
of a particular gene. Blue columns, experimental expression values (expression ratios in log2); red columns, simple moving averages of experimental log2 
values of 10 neighboring genes. Black arrows indicate position of ORFs analysed by qPCR as shown in Figure 6. 



chromosome using real time qPCR. The results (Figure 6) 
showed significant differences between gene copy num- 
bers in the left and right arms of chromosome XV in BR- 
RF indicating partial genome duplications. 

BR-F colonies up-regulate expression of genes located in 
HAST domains, which affect BR-F stress resistance 

Analysis of the differently expressed genes revealed that 
some are localized near telomeres. It has been hypothe- 
sized that transitions between different phenotypes 
are controlled by transcriptional silencing via specific 
chromatin-remodeling enzymes [18,19]. Thus, we com- 
pared the genes that were expressed differently in the BR-F, 
BR-S and BR-RF colonies with the list of genes expressed 
differently in a mutant containing a deletion of SIR2 [20], 



the gene encoding histone deacetylase. Histone deacetylase 
Sir2p has been shown to affect transcription at specific 
chromosomal locations near telomeres [21]. This analysis 
revealed a correlation between the genes down-regulated 
in the sir2A strain and the genes typically expressed in col- 
onies with a biofilm phenotype (29% of genes up-regulated 
in 4-day-old structured biofilm colonies). However, neither 
the structured colony morphology nor the domestication 
rate of the BR-F-5/r2 strain (in which both SIR2 alleles were 
deleted) was changed when compared with the BR-F strain 
(data not shown). This result implies that the function of 
Sir2p is not important for the formation of structured col- 
ony morphology. 

Next, we compared our set of genes with a database 
of genes specified to be under the control of another 



St'ovicek et al. BMC Genomics 2014, 15:136 
httpy/www.biomedcentral.com/l 471 -21 64/1 5/1 36 



Page 10 of 16 




2.32 ; 4.38 1.97 i 0.61 0.90 1.15 



Figure 6 Relative copy number of genes located on the left 
and right arm of chromosome XV in BR-RF versus BR-S. 

Comparison of relative copy number was performed by qPCR using 
cliromosomal DNA isolated from BR-RF and BR-S strains. Positions of 
examined genes on chromosome XV are marked by black arrows at 
Figure 5. Red, genes located at the left arm of the chromosome XV the 
expression of which is upregulated in BR-RF when compared with BR-S 
(and BR-F); blue, genes located at the right arm of the chromosome 
XV. Y axis, relative gene copy number (BR-RF versus BR-S); the exact values 
are stated below the ORFs. The mean of values from 3 independent 
experiments ± SD are shown; data significance (relative copy number of 
each left-arm-gene compared to relative copy number of each right-arm- 
gene) was determined using the two-tailed ttest. All P values varied 
between P < 0.05 and P < 0.005. 



deacetylase, Hdalp. Hdalp was shown to cause transcrip- 
tional repression in specifically localized subtelomeric 
chromosomal regions [22]. This analysis revealed a correl- 
ation between the genes localized to HAST domains, i.e., 
regions controlled by Hdalp histone deacetylase [22], and 
the genes up-regulated in BR-F colonies. In particular, 
18% of the genes up-regulated in BR-F versus BR-S col- 
onies, 19% of genes up-regulated in BR-F versus BR-RF 
colonies and 32% of genes up-regulated in BR-F versus 
both BR-RF and BR-S 4-day-old colonies were localized to 
these domains. In contrast, upregulation of HAST domain 
genes in BR-S and BR-RF colonies versus BR-F colonies 
was less than 6%. Thus, the upregulation of these genes in 
BR-F colonies indicated decreased activity of Hdalp com- 
pared to BR-S and BR-RF colonies. 

Genes in the HAST domains that were specifically up- 
regulated in the BR-F strain (Figure 7 A) included genes 
encoding various glucosidases (8 genes), proteins with 
mostly unknown function and predicted localization to 
the cell wall (11 genes) and various transporters (8 genes). 
These data indicated that the cells in the BR-F colonies 
could have specific cell wall and plasma membrane prop- 
erties compared to the BR-S and BR-RF colonies; however, 
these properties seemed unrelated to the structured 



biofilm colony phenotype. To test this prediction, we pre- 
pared the strain BR-S-hdal, with both alleles of the HDAl 
gene deleted, which should de-repress HAST domain 
genes. We then analyzed the morphology of colonies of 
this strain and the resistance to various stresses compared 
with the BR-F, BR-RF and BR-S strains (Figure 7B). The 
B^'S'hdal strain formed smooth colonies similar in 
morphology but slightly larger than colonies of the BR-S 
strain of the same age (Figure 7B). Similarly to the BR-F 
strain, the B^RS-hdal strain was significantly more resist- 
ant to the chitin binding dye calcofluor white, and it grew 
better at a higher temperature (37°C) than either the BR-S 
or the BR-RF strain. Thus, deletion of the HDAl gene in 
the BR-S strain increased its calcofluor white and 
temperature resistance to the same level as the BR-F strain, 
but it did not change the smooth colony morphology. 

Discussion 

Wild S, cerevisiae strains are able to effectively domesticate 
when transferred to favorable laboratory conditions [13,14]. 
Here, we show that the BR-S strain can switch back to a 
wild-type-like phenotype and start to form structured bio- 
film colonies following long-term exposure to starvation 
and stress conditions. In contrast to the original BR-F 
strain, which forms structured colonies with a uniform 
morphology, our isolated feral BR-S subclones form col- 
onies with different phenotypes. Thus, the increased pheno- 
typic variability of the feral BR-S subclones seems to be 
induced by environmental stress and starvation. These data 
agree with previous findings that non-isogenic S, cerevisiae 
strains isolated from various habitats can form differently 
structured colonies [14], similarly to pathogenic fungi [5]. 

The first feral subclones of the BR-S strain appear after 
approximately 30 days of slow growth in static cultures 
with limited nutrients. Later, the proportion of subclones 
increases up to 26% of the viable cells in some of the par- 
allel cultivations. At least between days 75 and 110, such 
an increase could be caused by an increased rate of 
phenotypic switching as well as by the preferential dying 
of smooth cells during the "loss of viability" phase of the 
population growth curve (Figure lA). The number of feral 
cells starts to decrease from day 115, i.e. 40 days later that 
in the case of domesticated cells. These findings suggest 
that feral subclones gain the ability to better survive star- 
vation than domesticated cells. Higher concentrations of a 
better carbon source (e.g., fructose) block the emergence 
of subclones forming structured colonies, thus indicating 
the importance of nutrient restriction for the evolution of 
feral subclones. Similarly, various stress factors promote 
phenotypic switching of C. albicans^ resulting in the pro- 
duction of cell variants that are more resistant to stresses 
[23]. Exposure of C. neoformans cells to various host 
defense mechanisms in vivo also leads to phenotypic 
modulation of an originally homogeneous population and 
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Figure 7 Derepression of HAST domain genes. A. Genes localized to HAST domains [22] up-regulated specifically in 4- and 7-day-old BR-F colonies 
compared to both BR-RF and BR-S colonies. Genes only expressed differently on the 4th day are in standard letters, genes expressed differently on both 
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to the subsequent prevalence of variants that are more re- 
sistant to the host immune response [6]. 

A BR-RF strain forms structured colonies with a 
morphology that is visually identical to that of the colonies 
formed by the original wild BR-F strain. The BR-RF strain, 
however, differs in several aspects from the BR-F strain. 
First, the BR-RF strain domesticates 4 times more effect- 
ively on fermentable media and 8 times more effectively 
on respiratory media than the BR-F strain. This finding 
suggests that the several months of starvation are suffi- 
cient to induce the formation of feral subclones but are 
not sufficient to make the feral subclones as stable as the 
original BR-F strain, which evolved over a much longer 
time in a natural setting. This finding also agrees with the 
results from some feral subclones that emerged during 
static cultivation (Figure IB); these subclones form less 
structured colony phenotypes that are even more unstable 
and quickly revert back to the domesticated phenotype. 



The data also demonstrate that the presence of a good fer- 
mentable carbon source such as glucose efficiently promotes 
the domestication of both strains (Figure IC) and prevents 
the switch to a structured phenotype as indicated above. 

Second, the colonies of the BR-RF strain differ from the 
colonies of both the BR-F and BR-S strains in certain 
physiological parameters, behaving as an "intermediate" be- 
tween the BR-F and BR-S colonies (Figure 2). Similarly to 
the BR-S strain, the BR-RF strain also forms oval cells and 
is not able to form the elongated cells and pseudohyphae 
that are typical of the BR-F strain. This finding supports 
the idea that the particular cell morphology is a dispensable 
factor in the formation of biofilm colony architecture [14]. 

On the other hand, both the formation of the ECM and 
the production of Flollp adhesin, the traits important for 
the formation of structured colony morphology [14], recover 
in BR-RF colonies. The length of the FLO 11 transcript, how- 
ever, differs from the original BR-F and the feral BR-RF 
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strains. Recombination of the central repeat domain of the 
yeast adhesin genes is frequently observed [24], which often 
leads to the generation of a protein with different and/or en- 
hanced adhesive function [25]. The change in the FLO 11 
gene of the BR-F strain most likely occurs during the do- 
mestication process; the FLO 11 transcript, although present 
at very low levels in BR-S colonies, is the same size in both 
the BR-S and BR-RF strains. Thus, it seems that rather than 
actual gene size and potentially protein size, the Flo lip ex- 
pression level is important in the context of structured col- 
ony formation. The altered FLO 11 transcript length also 
indicates that DNA recombination may occur frequently 
during phenotypic transitions of the BR-F strain. This con- 
clusion is further supported by the observed differences in 
the MW of the HMWGP protein component of the ECM 
(Figure 3), which could be caused by different posttransla- 
tional modifications or by changes in gene structure. In 
summary, these data show that the BR-RF strain, which was 
selected under the conditions of long-term nutrient 
deprivation and stress, develops some properties similar to 
the BR-F strain (e.g., the formation of structured biofilm col- 
ony morphology) but also maintains some properties of the 
BR-S strain (e.g., the formation of oval cells). 

Genome-wide transcriptome analysis allows us to deter- 
mine the typical gene expression profile for a structured 
colony phenotype compared with the smooth phenotype. 
Structured biofilm colonies (BR-F and BR-RF) both up- 
regulate and down-regulate distinct genes involved in cell 
wall remodeling as well as genes involved in the delivery 
and modification of secreted and cell wall/membrane-at- 
tached proteins. This finding indicates that, in addition to 
Flo lip adhesin, the carbohydrate composition of the cell 
wall, the appearance of specific surface proteins and pro- 
tein glycosylation may affect specific properties of the cells 
(such as the production of the ECM) and contribute to 
the formation of a biofilm colony. Ty transposon genes 
and several unknown genes that are expressed in biofilm 
colonies are located in subtelomeric regions. This suggests 
that transcription of genomic regions that are influenced 
by chromatin structure [26] differs in biofilm colonies 
compared to smooth colonies. In support, smooth colonies 
up-regulate the expression of a group of genes involved in 
histone modification and chromatin remodeling (Table 1). 

There is further evidence of the importance of chroma- 
tin remodeling in the diversification of properties of 
phenotypic variants derived from the BR-F strain. BR-F 
colonies (but not BR-RF colonies) express genes located in 
HAST domains, i.e., regions in the genome whose expres- 
sion is modulated by histone deacetylase Hdalp [22]. De- 
letion of the HDAl gene (and thus de-repression of 
HAST domain genes) significantiy increases the resistance 
of BR-S cells to both calcofluor white and high 
temperature, thereby reaching the level of resistance of 
the BR-F strain. Calcofluor white interferes with cell wall 



assembly [27], and strains with mutations that affect the 
cell wall composition are more sensitive to calcofluor 
white [28]. Thus, the resistance of the BR-F strain could 
be caused by expression of HAST domain genes encoding 
various cell wall and plasma membrane-related proteins. 
In support of persistent repression of HAST domain genes 
in BR-RF colonies, cells from BR-RF colonies behave simi- 
larly to BR-S cells. On the other hand, de-repression of 
HAST domain genes in BR-S-hdal colonies does not 
affect the smooth morphology of these colonies, indicating 
that HAST domain genes expressed in the BR-F strain are 
unrelated to the biofilm colony phenotype. However, other 
studies demonstrate the role of specific deacetylases, in- 
cluding Hdalp, in controlling the phenotypic transition 
of the "white-opaque" C. albicans colony [18,29]. 
Hence, epigenetic mechanisms mediated by chromatin 
remodeling enzymes have the capacity to regulate 
phenotypic heterogeneity. 

In addition to the expression characteristics specific to 
structured biofilm colonies, gene expression specific to 
BR-RF colonies is related to three regions of the BR-RF 
genome. This finding indicates that some chromosome re- 
arrangements occur during BR-RF strain formation; these 
could result in the duplication of certain parts of chromo- 
somes or in chromosomal rearrangement(s) relocating the 
regions to more actively transcribed locations of the gen- 
ome. Real time qPCR comparison of relative copy number 
of three genes located in the left arm of chromosome XV 
that upregulate expression in BR-RF {YOLllOW, 
YOL060C and YOL038W) with that of three genes located 
in the right arm of the same chromosome {YOR194C, 
YOR250C, YOR288C), indicates large duplications in the 
genome (Figure 6). This finding also agrees with previous 
observations that diverse conditions can evoke aneuploidy 
in yeast, which can directly impact gene expression, as de- 
tected at both the transcriptional and proteomic levels 
[30,31], and can also modulate the switch between differ- 
ent colony pheno types [32]. However, the observed 
chromosome site-specific expression characteristics of the 
BR-RF strain are unrelated to the structured colony 
phenotype, as only a few biofilm colony-specific genes are 
localized to these chromosomal regions. DNA re- 
arrangement events leading to karyotypic instability are 
also observed during phenotypic transitions of certain 
C. albicans [33] and C. neoformans strains [8]. These 
previously observed changes are not correlated with 
phenotypic variability and are irreversible. Thus, it 
seems that DNA rearrangements resulting in particular 
gene expression changes are not the cause of the BR-F/ 
BR-S/BR-RF phenotypic switch, which is in contrast to 
recent findings in other wild strains [32]. Instead, these 
rearrangements in the BR-RF strain could be a side- 
product of the increased rate of recombination events 
that occur during the phenotype transitions. 
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Conclusions 

We describe the phenotypic, genomic and gene expression 
differences among three S. cerevisiae strains (wild, domesti- 
cated and feral) that emerged by phenotypic switching 
under diverse environmental conditions. These strains form 
colonies with varying complexity and gain diverse features 
including altered resistance to stress. We show that the wide 
variability of natural S, cerevisiae strains can be further po- 
tentiated under stressful environmental conditions. We 
identify genes specifically expressed in the structured colony 
phenotype, such as genes that specifically affect the compos- 
ition of cell surface structures; these genes provide clues as 
to the specific processes involved in the formation of struc- 
tured colony architecture. We also document a role of 
Hdalp histone deacetylase in strain resistance to stress. In 
addition, our data indicate that genomic rearrangement(s) 
occur in the feral strain that are unrelated to the phenotypic 
switch. These changes contribute to expression characteris- 
tics of the feral strain that are distinct from both its ances- 
tors, the wild and domesticated strains. 

In summary, our findings show that extensive phenotypic 
modulation occurs in wild S, cerevisiae strains. Phenotypic 
switching can be controlled by various mechanisms and 
enables flexible adaptation of these strains to a particular 
environment in the wild. The comparison of the three 
originally isogenic strains BR-F, BR-S and BR-RF re- 
vealed many significant gene expression and phenotypic 
differences, but only a few of them seem to be related 
to the formation of structured biofilm colonies. In 
addition, the differences observed between BR-F and 
BR-RF showed that reverted structured colony morph- 
ology does not necessarily mean reversion of other 
strain properties. In other words, strains with different 
properties such as BR-F and BR-RF can form structured 
colonies. Thus, an important question remains as to 
whether there is a single dominant mechanism that 
forces yeast strains to form structured colonies inde- 
pendently of other properties or, more probably, a var- 
iety of mechanisms that support the formation of 
structured colony morphology. 

Methods 

Strains and media 

All strains and their derivatives are listed in Additional file 
2: Table S2. The wild S, cerevisiae strain BR-F, originally 
collected from the lake in south Slovakia [34] is from the 
collection of the Chemical Institute of the Academy of Sci- 
ence, Bratislava, Slovak Republic (collection number CCY 
21-4-97). The BR-RF and BR-S strains are isogenic to the 
BR-F strain and arose from phenotypic switching. Col- 
onies were grown on GMA (3% glycerol, 1% yeast extract, 
2% agar) or YEGA (2% glucose, 1% yeast extract, 2% agar) 
media at 28°C unless otherwise indicated. For static culti- 
vations, we used liquid minimal medium (MM) (0.5% 



(NH4)2S04, 0.1% KH2PO4, 0.05% MgS04, pH 5-6.5) sup- 
plemented with ethanol (0.1%, 0.5%, 1%, 2% and 3% v/v), 
fructose (0.2%, 0.5%, 1% and 2%) or NaNOs (0.01% w/v). 
For selective plates, media were supplemented with 
200 mg/1 G418 or 100 mg/1 nourseothricin. 

Static cultivations, determination of CFU and percentage 
of cells that switched to form biofilm colonies 

The BR-S strain was inoculated into MM medium at a 
concentration of approximately 5x10^ cells per ml. The cell 
culture was split into 10 ml portions in glass tubes (3 to 10 
parallels for each particular concentration of the supple- 
ments) and statically cultivated at 28°C. During long-term 
cultivation, the culture was mixed at particular time- 
points. After appropriate dilution, 100 \A aliquots were 
plated on GMA to determine the number of CFU and the 
proportion of cells forming structured biofilm colonies. 

Construction of the strains 

Flollp-GFP strains containing the FLO 11 gene fused with 
the GFP gene were constructed via a previously described 
procedure [17]. The BR-RF strain with an artificial promoter 
Pgali-GFP cassette integrated in the genome instead of the 
HISS gene was constructed as previously described [16]. 
Gene knock-outs were performed by transforming the cells 
with deletion cassettes generated by PCR using the primers 
and plasmids listed in Additional file 2: Table S3. Yeast cells 
were transformed according to a published protocol [35]. 

Determination of colony biomass 

The biomass of colonies was estimated as their wet 
weight. To determine dry biomass, the wet biomass was 
dried in a Speed- Vac for at least 6 h. The water content 
was calculated as the difference between the wet and dry 
biomasses. The data were presented as the average of 
three independent experiments ± SD (standard deviation). 

Determination of switching frequency 

The occurrence of switching from the BR-F and BR-RF 
to the BR-S colony phenotype was determined by visu- 
ally scoring the colonies of different ages growing either 
on GMA or on YEGA media. To calculate the switching 
frequencies, 100-200 BR-F and BR-RF cells were plated 
on GMA plates, and at least 5000 colonies were scored. 
The percentage of colonies of the total CFU with 
smooth morphology was determined. 

Fluorescence microscopy of cells and colony imaging 

Cells were examined under a Leica DMR microscope 
equipped with a lOOx/1.3 oil objective and a GFP filter or 
Nomarski contrast and photographed with a ProgRes® 
jyjpcooi QQY>^ camera (Jenoptik, Germany). Colony images 
were captured in incident light. A ProgRes® CT3 CMOS 
camera with a Navitar objective, Fiber-Lite PL-800 
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illumination system and NIS Elements software (Labora- 
tory Imaging) were used. 

Sensitivity of the strains to calcofluor white 

The strain sensitivity was tested by the drop assay. 
Four-day-old colonies were harvested to make a sus- 
pension containing 2x10^ cells per ml. Next, 5 \i\ 
drops of 10-fold serial dilutions were applied to GMA 
plates supplemented with 10 mg/1 or 200 mg/1 calco- 
fluor white. Plates were scored after 3 days of incuba- 
tion at 28°C or 37°C. 

Northern blot analysis 

Total RNA from the cells of colonies grown on GMA 
plates was isolated by the hot phenol method as previ- 
ously described [14]. For northern blots, 12 [ig of total 
RNA was separated on a 1.5% agarose gel, transferred to a 
positively charged nitrocellulose membrane (Hybond-XL, 
Amersham Bioscience) and exposed to a labeled probe. 
The radioactive signal was visualized on Fuji X-ray film. 
The DNA probe for the RDN18 gene was a complete ORF 
of the gene prepared by PGR reaction. For the FLO 11 
probe, a PGR fragment corresponding to the last 1382 bp 
of the FLO 11 gene was used. The [a-^^P] dGTP-labeled 
probes were obtained by random priming using the Deca- 
Prime II Kit (Ambion). 

MIcroarray analysis 

The mRNA was isolated from total RNA using a Micro-Fast 
track 2.0 Kit (Invitrogen). Reverse transcription was per- 
formed from 2-4 (ig of mRNA as previously described [36]. 
Labeled cDNA was applied to Yeast 6.4 K Array (Y6.4 K) 
double spotted ORF slides (University Health network, 
Toronto, Ganada) following the manufacturer s instruc- 
tions and incubated overnight at 37°G. The arrays were 
washed according to the manufacturers protocol and 
scanned using an Axon fluorescent scanner and the 
GenePix software. Three double-genome microarrays 
(biological replicates) were used. Spots were detected 
using the TIGR Spotfinder (TM4.org) software and normal- 
ized in TIGR MIDAS (TM4.org) using the LOWESS 
method. Differentially expressed genes were selected using a 
combined fold-change and ANOVA analysis using MeV 
software (TM4.org). Genes exhibiting an at least 1.5-fold 
change (log2 > 0.585) in their average expression and simul- 
taneously a p-value < 0.05 were considered to be differen- 
tially expressed. 

Extracellular material extraction, SDS-PAGE and detection 
of protein glycosylation 

Extracellular material was extracted as previously described 
[14]. After SDS-PAGE separation, extracellular proteins 
were stained with silver [37] or transferred to a PVDF 



membrane. Glycoproteins were visualized by sequential 
incubation in concanavalin A, horseradish peroxidase 
and chloronaphthol H2O2 solution [38]. 

Two-photon excitation confocal microscopy (2P-CM) 

Vertical transverse cross-sections of colonies were prepared 
and their side-views obtained by 2P-CM as described 
previously [39]. The cross-sections were stained with 
Concanavalin A conjugated with Alexa Fluor 488 
(ConA-AF; 30 (ig/ml). Alternatively, GFP fluorescence 
was monitored. Images were acquired with a Leica TGS 
SP2 AOBS MP confocal laser scanning microscope fit- 
ted with a mode-locked Ti: Sapphire Chameleon Ultra 
laser (Coherent Inc., Santa Clara, CA, USA) for two- 
photon excitation as previously described [16]. 

Real time quantitative PGR (qPCR) analysis of genomic DNA 

Three target genes were selected from those in the left arm 
of chromosome XV with apparently upregulated expression 
in strain BR-RF compared with BR-S and three reference 
genes were selected from those in the right arm of chromo- 
some XV with no significant difference in expression 
(Figure 5). Primers (Additional file 2: Table S4) were de- 
signed using the GenScript real time qPCR primer design 
tool (https://www.genscript.com/ssl-bin/app/primer). Stipu- 
lated parameters included an amplicon length of 90-110 bp 
and melting temperature of 60-65°C. Primers were selected 
with a GC content of 40-60%, a length of 19-21 bp 
and minimal self-complementarity when analyzed using 
OligoCalc (http://www.basic.northwestern.edu/biotools/oli- 
gocalc.html [40]). DNA extraction was carried out according 
to the method of [41]. Real time qPCR was carried out on 
the light cycler 480 II (Roche, Basel) using FAST START 
Sybr Green (Roche, Basel), forward and reverse primers 
(final concentration 500 nM) and genomic DNA (20 ng per 
20 (il reaction) from strain BR-RF or strain BR-S (or water 
in negative controls). An initial melting/activation step at 
95°C for 5 min was followed by 45 cycles of melting, anneal- 
ing and amplification (95°C for 30 s, 56°C for 30 s and 72°C 
for 30 s) with a fluorescence measurement (465-510 nm) at 
the end of each amplification step. Raw data were analyzed 
using the comparative delta CT method [42]. The relative 
gene copy number is 2^^^, where ACT = (CTbr_s-CTbr-rf). 

In silico analysis of gene sets obtained by microarray 
comparison 

Genes were classified into functional categories using the 
particular gene description in the Saccharomyces genome 
database available at http://www.yeastgenome.org/. En- 
richment of the gene sets at particular chromosomal re- 
gions was determined by Positional Gene Enrichment 
analysis using the online tool available at http: //homes. 
esat.kuleuven.be/~bioiuser/pge [43] . 
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Availability of supporting data 

The data sets supporting the results of this article are 
available in the NCBI GEO repository (accession num- 
ber GSE40625, http://www.ncbi.nlm.nih.gov/geo/query/ 
acc.cgi?acc=gse40625). 

Additional files 



Additional file 1: Table 51. List of genes differentially expressed 
comparing transcriptomes of cells from 4- and 7- day-old BR-F, BR-RF and 
BR-S colonies. Differentially expressed genes were selected on the basis 
of combined criteria of fold-change (0.585 on log2 scale) and p-value 
(<0.05). Gene sets obtained from comparison of each phenotype couple 
(BR-F X BR-S, BR-F x BR-RF and BR-S x BR-RF) for both 4- and 7-day-old 
colonies are on a separate xls list. 

Additional file 2: Table 52. List of strains used in this study. Table S3. 
List of primers and plasmids used in the strain construction. Table S4. 
List of primers used in qPCR. 
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